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I , We study the chiral properties and the vahdity of perturbation theory for 

\ domain wah fermions in quenched lattice QCD at (3 = 6.0. The exphcit chiral 

symmetry breaking term in the axial Ward-Takahashi identity is found to be 

■ very small already at Ng = 10, where A^^ is the size of the fifth dimension, 

■ and its behavior seems consistent with an exponential decay in Ng within the 
^) , limited range of Ng we explore. From the fact that the critical quark mass, 

\ at which the pion mass vanishes as in the case of the ordinary Wilson-type 

I fermion, exists at finite Ng, we point out that this may be a signal of the 

parity broken phase and investigate the possible existence of such a phase in 
this model at finite Ng. The p and vr meson decay constants obtained from 
, the four-dimensional local currents with the one-loop renormalization factor 

^ \ show a good agreement with those obtained from the conserved currents. 
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I. INTRODUCTION 



Domain wall fermion model is a 5-dimensional Wilson fermion with free boundaries 
in the fifth dimension. At each of the boundaries either of the left- and right-handed chiral 
mode is expected to be localized exponentially, and hence the mixing between them, which 
represents the explicit chiral symmetry breaking effects, is suppressed exponentially. In the 
large Ng limit, where Ng is the size of the fifth direction, domain wall QCD (DWQCD) 
would have the desirable features, such as (i) no fine tuning for the chiral limit, (ii) O(a^) 
scaling violation and (iii) no mixing between three- and four-quark operators with different 
chiralities. It is shown that these expectations are realized within the perturbation theory 
up to the one- loop level in the limit Ng ^ oo Pioneering numerical studies also seem 

to support these expectations 

In this paper we focus on two issues. One is the chiral properties of DWQCD at the finite 
Ng. It is of great importance to understand correctly how the ideal chiral properties are 
realized nonperturbatively as A^^ increases. We investigate the Ng dependence of the explicit 
chiral symmetry breaking term in the axial Ward-Takahashi identity, which are expected 
to be exponentially suppressed as the A^^ increases. Correspondingly, at finite A^s, the pion 
remains massive in the chiral limit. We notice, however, that the pion mass at finite A^s may 
vanish at some non-zero value of quark mass, which we call the critical quark mass, in a 
similar manner that the pion mass in the Wilson fermion formulation vanishes at the critical 
hopping parameter. In the case of the Wilson fermion formulation, it is well known that 



the massless pion implies the existence of a parity(-fiavor) broken phase [|TT|. The analysis 
of two-dimensional Gross-Neveu model with domain-wall fermions indeed demonstrates the 
existence of the parity broken phase at the finite A^^, and the negative critical quark mass 
goes to zero as Ng increases |jl^,[l3|. In this paper we explore the negative quark mass 
region to examine the possible existence of the parity broken phase in the four- dimensional 
DWQCD. 

Another important issue in this paper is the validity of the perturbation theory for 
DWQCD. Recently we calculated the renormalization factors for the bilinear operators 
and three- and four-quark operators [0], for the latter of which we find the multiplicative 
renormalizability at infinite A^^ as opposed to the Wilson fermion case. We test the validity 
of the perturbation theory comparing the p and tt meson decay constants extracted from 
the local currents with those from the conserved ones. 

This paper is organized as follows. In Sec. |T| we present the details of simulation including 
the action of domain wall fermion. The chiral properties of DWQCD are investigated in 



Sec. |T|, where we also examine the existence of the parity broken phase. In Sec. |3 we test 
the validity of the perturbation theory for DWQCD using the meson decay constants. Our 
conclusions and discussions are summarized in Sec. 0. 

II. DETAILS OF NUMERICAL SIMULATION 
A. Action 

The domain wall fermion action is written as 
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Sf = - ^{^^s)Ddwf{x,s]y,s')i>{y,s') , (1) 

Ddu,f{x, s; y, s') = (5 - M)S^^ySs,s' - D^{x, y)5s,s' - D^{s, s')S^^y , (2) 
1 
2 



D\x, y) = l-\{l- -if,)U^J^+^^y + (1 + lMJx-[,,y] , (3) 



Pl52,s' {s = 1) 

D\s, s') = { Pl6s+i,s' +PrSs-i,s' {l<s< N,) (4) 

PrSns-i,s' {s = Ns) 

with domain wall height M and the extent of the fifth dimension Ng. The right- and left- 
handed projection operators are defined by Pr^l = (1 ± 75)/2- The chiral zero mode is 
supposed to appear for M > 0. The index x, y stands for four-dimensional lattice sites, 
while the index s, s' for the fifth direction runs from 1 to iV^. 

The quark fields on the four-dimensional space-time are given by combinations of the 
fermion fields at the boundaries, 

q{x) = PL^Ijix,l) + PR^Ijix,Ns) , (5) 
q{x) = tPix, 1)Pr + tPix, Ns)Pl , (6) 

where the left(right) chiral zero mode is well localized at s = l(A^s) for < M < 2 in the 
free theory. 

The explicit bare quark mass is introduced as 

Sf^Sf-Y. mfq{x)q{x) , (7) 

X 

which provides the quarks in eqs.(^ and with a Dirac mass rhf. For the free quark case, 
we obtain frif = M{2 — M)mf in the Ng — > oo limit 0. 

The action (|I|) is invariant under global SU(A^/) symmetry, which yields the five- 
dimensional conserved current 0, 



ip{x, s)(-l + -i^,)U^^^—i){x + /i, s) + ij{x + /i, s)(l + l^,)Ul^^—i){x, s) 



(8) 



where A'^/2 (a = 1, . . . , N'j — 1) are generators of SU(A^/) group. In terms of this current we 
can define the conserved vector and axial currents, 

Ns 

vr''i^) = Y.j>^^)^ (9) 

s=l 

^a,con^^) = - E Sign ( Y - " + 2 ' 



s=l 



We should note that A'J^^"'^ is not exactly conserved in the finite Ns contrary to the vector 
case: the divergence of satisfies 

V-^A';:^°\x)0{y)) = 2mf Ux)^,^q{x)0{y)) + 2 lr,,{x)0{y)) (11) 



with 

Jt,{x) = -ij{x, NJ2)Pl^^{x, Ns/2 + 1) + tp{x, Ns/2 + 1)pJl.^{x, iV,/2), (12) 

where is the backward difference operator: V~f{x) = f{x) — f{x — fi) and O is some 
operator. We will set x y in our numerical simulation. The contribution of J^q{x) 
in eq.(0) is expected to be exponentially suppressed as the A^^ increases for the physical 
operator Ophys.(^, 



B. Simulation parameters 

The simulation is carried out with the plaquette gauge action aX (3 = 6.0 on 16^ x 32 x A^^ 
lattices in quenched QCD. Table | summarizes our run parameters. We generated gauge 
configurations on a 16^ x 32 lattice with the five-hit pseudo heat-bath algorithm. Each 
configuration is separated by 2000 sweeps after 20000 sweeps for the initial thermalization. 
In DWQCD the same gauge configuration is assigned on each layer of the fifth direction. 
The p meson mass at the chiral limit (m/ = 0) is used to determine the lattice spacing, 
which gives a^^ ~ 2GeV with rup = 770MeV as input while slightly depending on M and 
Ng (see Sec. PTTI ). The mean- field estimate for the optimal value of M gives 0] 

M= 1 + 4(1 -m) = 1.819 , (13) 

employing u = l/{8Kc) with = 0.1572 for the Wilson fermion at /3 = 6.0. 

The quark propagator is obtained by solving the inverse of the quark matrix (0) with 
the unit wall source without gauge fixing, 

E Y: Dtfix, s; y, s')G'^{y, s') = B\x, s) , (14) 

y,s' b 

B'^ix, s) = ^Y. ^acSt,o E [^s,iPr + Ss,nM , (15) 

'iVs c S' 

where x = {x, t) with spatial coordinates x and time coordinate t; Vs is the spatial volume 
and a, b, c denote the color indices. The four- dimensional quark propagator is constructed 
with a linear combination of G^, 

G,ix) = PlG^x, 1) + PrG^x, Ns) . (16) 

To solve the linear equation ([l^) , we employ the conjugate gradient method with even/odd 
preconditioning In the case of nif = 0.025, about 500 iterations are needed to satisfy 



a convergence criterion which requires Irp/IG^p < 10~^^ with the residual r. Neither 
minimal residue(MR) nor BiCGStab converges within a practical number of iterations for 
gauge configurations at /3 ~ 6.0. Throughout this paper, statistical errors are estimated by 
the single elimination jackknife method. 
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III. CHIRAL PROPERTIES OF DWQCD 



A. Restoration of chiral symmetry 



We investigate the existence of the chiral zero mode expected in the hmit of mj ^ and 
Ng oo using the pion mass at the chiral limit [rrif = 0) and the anomalous contribution 
Jsq in the PC AC relation (0). The pion mass mT^{mf, Ng, M) is extracted from a global fit 
of the two-point function 

J2{Mx,t)Jl{y,0)) (17) 

x,y 

employing the fitting function 

where Jt,{x) = q{x)'^^q{x). We choose the fitting range to be 9 < t < 23 after taking account 
of the time reversal symmetry t ^ 32 — t. 

Figure |1| shows the m^a dependence of m^. We observe a good linear behavior for the pion 
mass squared as a function oi nif. Employing a linear extrapolation, we obtain the value of 

at = 0, which is shown in Fig. |^ as a function of M. We find that the minimum point 
of the pion mass is around M = 1.819, which shows the mean- field estimation of optimal 
M in (|T3p is reasonably accurate. Furthermore, extending the linear extrapolation into the 
negative rrif region, as is done for the ordinary Wilson- type fermion, we can also extract the 
critical quark mass mdNs, M) at which the pion mass squared vanishes. As we expected, 
both at = and mc{Ns, M) decreases monotonically as Ng increases. In Fig. ^ we 
plot ml for nifa =0.075, 0.050, 0.025 and as a function of Ns, where m^imf = 0) is the 
extrapolated value. The magnitude of at m/ = seems to diminish exponentially in 
Ng. In order to confirm that at mj = indeed vanishes exponentially in Ng, one has to 
further increase A^^. However, we do not attempt such an investigation at the present lattice 
size, since m^(mj = 0) at Ng = 10 is already as small as that for the Nambu-Goldstone 
pion of the Kogut-Susskind fermion, which gives an estimate of the finite size effect to the 
would-be massless pion at the same P and the spatial lattice size. 

It is reported, at /? = 6.0 on a 16^ x 32 lattice, that = 0.014(2) in the chiral limit 
for Ng = 16 and M = 1.8 [|1^, which is roughly equal to our value of = 0.0093(54) for 
Ng = 10 and M = 1.819. 

The anomalous contribution due to J^q in eq. ( pT]) is evaluated by calculating the ratio of 
two-point functions. 



E.^,,^(j.(x,t)4(y,0)) (OlJ.lvr)' 

where we employ a constant fit over the range 8 < t < 24. This quantity measures contri- 
bution of the explicit chiral symmetry breaking effects in the PCAC relation: 



0|V-A?°"|7r 



(0|2J5g|7r) 

(0|J,|7r) (0|J,|7r) 



- 2m/. (20) 
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In Fig. 1^ we plot the results for (0| J5g|7r)/(0| JTrlvr) as a function of m/a. The data shows 
little nif dependence at each Ng. We obtain the value of (0| J5q|7r)/(0| J^|7r) at rrif = 
by extrapolating the data linearly to m/ = 0. Figure |^ illustrates the A^^ dependence of 
(0| J5g|7r)/(0| J^Itt) at nifa =0.075, 0.050, 0.025, 0, where we also plot the results for \mc\ 
for comparison. We find that the magnitude of (0| J5g|7r)/(0| Jn-|7r) decreases exponentially 
as Ng increases. The same situation is observed for |mc|. These features suggest that the 
chiral symmetry would be restored in the limit of m/ — >■ and Ng —>■ oo at (3 = 6.0. 

Note that \mc\ and (0| Js^lvr)/ (0| J^-Itt) are expected to be consistent, if the latter has no 
rrif dependence as seen in Fig. ^. However we observe a slight deviation between them in 
Fig. 1^. We point out two possibilities to yield this inconsistency. One is the finite size effects 
on the pion mass and the other is the quenched chiral logarithm plaguing the pion mass in 
the quenched approximation. We should note that the quantity (0| J5q|7r)/(0| JttItt) is free 
from both systematic contaminations, which assures that (0| Js^l vr)/ (017,^1 vr) is superior to 
\mc\ to measure the remnant of the chiral symmetry breaking effects in DWQCD. 

We briefly mention the determination of the lattice spacing a from the p meson mass. 
Figure^ shows the A^^ dependence of m^a at = (circles) and at = rric (squares) with 
M = 1.819, which are obtained by a linear extrapolation from the data. For small Ns{= 4), 
rripa at my = differs from the one at mj = rric beyond statistical errors, while the difference 
almost vanishes at A^^ = 10. This is true for other values of M: within statistical errors the 
values of rripa at ruf = and mj = rric are consistent with each other for M = 1.7 — 1.9 
at Ng = 10. Furthermore the M dependence of rup itself is mild over this range of M at 
Ng = 10. /^From the value of rripa at rrif = with M = 1.819 and Ng = 10, we estimate 
a-^ = 2.02(17)GeV from nip = 770MeV as an input. 



B. Parity broken phase 

In the Wilson fermion formalism, the critical quark mass is identified with the second 
order phase transition point between the (spontaneously) parity broken phase and the parity 
symmetric phase [0. The massless pion is understood as the massless particle associated 
with the second order phase transition. Recently, the study of the two-dimensional Gross- 
Neveu model with domain wall fermions pointed out that this picture is true even for this 
model at finite 



T^,|T3| . In this section we examine whether the massless pion at m/ 



rrir. 



at finite A^^ can be understood in this way. 

For the finite Ng, the parity broken phase may exist in negative rnj regions. From the 
study of the two-dimensional Gross-Neveu model with domain wall fermions |T3] , we expect 



the phase diagram in Fig. |^ for four-dimensional DWQCD. Each of five critical points, 
where the = line is touched on g'^ = line atm/ = -(l-M + 2/)^" (/ = 0, 1, 2, 3, 4), 
corresponds to the massless particle pole of momentum 



p^(/ = 0) = (0,0,0,0), 



TT 



7r 



p,{l = 1) = (-, 0, 0, 0), (0, -, 0, 0), (0, 0, -, 0), (0, 0, 0, - 
a a a a 

,7r TT , ,vr TT , ,7r vr. 
p^il = 2) = (-, -, 0, 0), (-, 0, -, 0), (-, 0, 0, -) 
a a a a a a 

7T TT . , vr vr , , vr vr , 
,(0,-,-,0),(0,-,0,-),(0,0,-,-) 
a a a a a a 



(21) 
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,7r TT TT . ,7r TT TT. .TT TT TT . , TT TT TT . 

p^il = 3) = (-, -, -, 0), (-, -, 0, -), (-, 0, -, -), (0, -, -, -), 
a a a a a a a a a a a a 

/, .\ ,TT n IT TT 

Pm(^ = 4) = (-,-,-,-)• 

a a a a 

Since the critical point for the p^{l = 0) mode converges to = for < M < 2 as 
Ns increases, while other four critical points move rapidly to |m/| = oo, we can expect 
one small < 1) critical point and four large > 1) critical points at finite Ng in 

four-dimensional DWQCD. 

Based on these speculation, let us investigate existence of the "cusp" (branch) correspond- 
ing the Pfi{l = 0) mode, which provides an evidence of the parity broken phase. We calculate 
pion masses at m/a = -0.120, -0.100, -0.093, -0.086, -0.080 for A^", = 4 and M = 1.819. 
Without adding the external field which couples to the parity broken order parameter, one 
can not perform the correct simulation in the parity broken phase. To avoid the direct sim- 
ulation in the parity broken phase, we take rather large magnitudes for negative m/ values 
so that the system goes into another symmetric phase, and try to find another side of the 
phase boundary (critical quark mass) from this phase. In Fig. ^ we plotted the tt meson 
mass squared function of m/. Extrapolations of ml to zero both from positive 

and negative ruf, where four largest ruf are used for negative ruf, indicate that a parity 
broken phase may exist around mja ^ —0.03. We can interpret that Fig. represents the 
portion in Fig. |^ denoted by dashed segment. It is noted that the large errors for negative 
rrif are caused by the fact that the pion propagators form peculiar shapes similar to "W" 
character, which has been often observed near the another side of the parity broken phase 
for the Wilson fermion |lH . 



IV. MESON DECAY CONSTANTS 

A. Perturbative renormalization factors 

The p meson decay constant fp is defined by 

{0\V,\p)=e,rnl/fp, (22) 

where nip is the p meson mass and is the polarization vector. We have two choices for 
the vector current V^: the exactly conserved current in eq.(^ and the four- dimensional local 
current V^(a;) = q{x)'jpq{x) multiphed by the renormalization factor Zy. 
In a similar way, the n meson decay constant f-,, is given by 

{0\Ap\7T) = ppU , (23) 

where Pp denotes the momentum of the pion and Ap is either the (almost) conserved current 
in eq.(|TU|) or the four-dimensional local current AJ^^x) = q{x)'^p'y^q{x) multiplied by the 
renormalization factor Z^. 

Perturbative technique in DWQCD is already applied to evaluation of the renormaliza- 
tion factors for the bilinear operators and three- and four-quark operators consisting 
of four-dimensional quarks in eqs.(^ and (j|). It is shown that Zs = Zp = Z^ and 
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Zy = Za at A^s — > oo for the bilinear operators, which is expected in the case that the chi- 
ral Ward-Takahashi identities hold exactly [H- Another desirable feature is that the three- 
and four-quark operators at Ng — > oo can be renormalized without mixing between op- 
erators with different chiralities, as opposed to the Wilson fermion case 0. The peculiar 
feature in the renormalization of DWQCD, however, is an appearance of the overlap factor 
(1 — |1 — M\'^)Zy^ for the four- dimensional quark fields. The one-loop coefficient of 
in is of 0(10 — 10^) for |1 — M|^0.1 without mean-field improvement, which could 
be problematic because the present simulation is done for M > 1.7. It is shown in Ref. 
0, however, that the magnitude of the one-loop coefficient for Z^j is reduced to 0(1) for 
|1 — M|^0.1 after the mean field replacement from M to M = M + 4(u — 1). Therefore, com- 
paring the p{7i) meson decay constant fp{fn) obtained from ZyV^(ZAA^ = ZvAJJ with that 
from V^°'^{A^°^) gives a good testing ground for the validity of the (mean field improved) 
perturbation theory. 

Following the notation of Ref. 0, the renormalization factor of is written as 

where 

wo{M) = 1 - M, (25) 
Z^M) = 1 + -^Cpz^iM), (26) 

Zv{M) = 1 + J^CpZviM), (27) 

without mean-field improvement. Cp = {N"^ — 1)/{2N) is the quadratic Casimir in SU(A^) 
gauge group. In the case of M = 1.819, we find 1 - = 0.32924, z^{M) = 250.36 and 
zv{M) = —17.760 using the results in Ref. 0. The large corrections of 1/(1 — Wq) and 
Zw could jeopardize the efficiency of the perturbation theory. However, once a mean-field 
improvement is employed, the renormalization factors are reexpressed as 

V™" = J ^uZ^^(M)V;^ (28) 

(1 - ^o(M)2)ZMF(M) ^ ^ ' ^ ^ 

with 

M = M + 4(m-1), (29) 

wo{M) = 1 - M, (30) 

^r(M) = l + Yf^C^F^r(M), (31) 

Zr(M) = 1 + r|^O^^r(M). (32) 

With the use of m = 0.87781 = v^, where P is the expectation value of the plaquette, we 
obtain M = 1.3302, 1 - Wo{My = 0.89095, z^^{M) = 7.327 and z^^{M) = -7.3744 from 
the results in Ref. It is remarkable that the perturbative corrections are drastically 
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reduced. We determine the coupling constant at the scale 1/a in the MS scheme with the 
aid of the mean- field improvement, 

^ = - 0.13486. (33) 



Incorporating g^{l/a) in eqs. (pi]) and (P^, we finally obtain the value of the renormaliza- 
tion factor with the mean- field improvement, 

^ -^uZ^^{M) = 0.7504. (34) 



w;o(M)2)zr(M) 



If we employ u = l/{8Kc) with Kc = 0.1572 for the mean field improvement, we obtain 
0.7147 instead of 0.7504 in eg. (|3^ . Since their difference is smaller than the statistical errors 
of the p and vr meson decay constants, we use the value of eq. (0) for the renormalization 
factor in the next subsection. 



B. p and vr meson decay constants 

We calculate fp from the ratio of two point functions, 

E.-.-(jp(x,t)j;(y,o) 



{0\JM {0\Vr'\p) ^ (35) 

where Jp{x) = q{x)'yiq{x) {i = 1,2,3) is the interpolating field for the p meson. The 
amplitude {0\Jp\p) is obtained by fitting the two-point function with 

E {Mx, t)Jl{y, 0)) = y j^^^y (e-™"* + e-''^^^-*)) , (36) 

x,y P 

where the fitting range is chosen to be 8 < t < 24. 

In Fig.^ we plot the p meson decay constants obtained from V^^ and V^°'^ as a function of 
ntf. We observe that both results are consistent once the perturbative corrections are applied 
to V^. This is an encouraging evidence to show the validity of the perturbation theory with 
the mean-field improvement. Still, the scaling violation effects should be checked. The value 
at the chiral limit = 0, which is obtained by extrapolating the results for the conserved 
current linearly, is found to be slightly smaller than the experimental value. 

Let us turn to the vr meson decay constant fj^. We obtain /^r from correlation functions of 
the axial vector current Ap{x) and the interpolating field for the vr meson J^(x) = q{x)'^^q{x), 



E (X(x, t)Y\y, 0)) = Vs ' (e-"^-* + e-^^^^^-t)^ ^ (37) 



x,y 



where X^Y = J^, ^4. The fitting range is chosen to be 9 < t < 23. 

Figure ^ shows a comparison of the tt meson decay constants obtained from and Af^. 



We find the same situation as in the case of the p meson decay constant: the perturbative 
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corrections sufficiently compensate the difference between Al and A'^'^. Linear extrapolation 
of the results for the conserved current gives the value at the chiral limit mj = 0, which is 
roughly consistent with the experimental value. We also show the M dependence of at 
ruf = in Fig. |ll|. Although the results slightly depend on the choice of M, the differences 
are smaller than current statistical errors. 



V. CONCLUSIONS AND DISCUSSIONS 

We have investigated the chiral properties of DWQCD by measuring the pion mass and 
the explicit chiral symmetry breaking term in the PCAC relation. Their Ng dependence up 
to Ng = 10 seem to be consistent with exponential decay in Ng, which indicates that the 
chiral symmetry is restored at A^^ cxd limit, as expected. More extensive study with larger 
Ng, however, is required to con&m this conclusion. 

We also calculate the pion mass at negative to examine whether there exists the 
parity broken phase in DWQCD. Our result is consistent with the existence of the parity 
broken phase. 

The validity of the perturbation theory in DWQCD is tested by calculating the p and vr 
meson decay constants from the four- dimensional local currents and the conserved ones. 
We find that the difference between both currents is sufficiently compensated with the 
perturbative renormalization factor up to the one-loop level with the aid of the mean-field 
improvement. This is an encouraging result, though the magnitude of the scaling violation 
on these quantities should be checked further. 
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FIG. 1. Pion mass squared as a function of m/a at M = 1.819 and Ng = 4,6, 10. Solid lines 
show linear fits. 
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FIG. 2. Pion mass squared extrapolated to m/ — > as a function of M at A^s = 4, 6, 10. 
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FIG. 3. Pion mass squared as a function of Ng &t M = 1.819. Filled square shows the value 
for the Nambu-Goldstone pion of the Kogut-Susskind fermion at /? = 6.0 on a 16^ x 40 lattice. 
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FIG. 4. (0| J5q|7r)/(0| J7r|vr) extrapolated as a function of m/a at M = 1.819 with Ng = 4,6,1 
Solid lines show linear fits. 
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FIG. 6. p meson mass extrapolated to = and m/ = mc as a function of Ng at M = 1.819 . 



14 




FIG. 7. Schematic phase diagram in {g^,mf) plane for A^^ =even. See text for dashed segment. 
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FIG. 8. in lattice unit as a function of m/a at M = 1.819 and Ng = 4. Solid lines show 
linear extrapolations. 
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FIG. 9. fp as a function of m/a at M = 1.819 and Ng = 10 together with the experimental 
value. Data points for ZyV^ (filled squares) are slightly displaced horizontally for clarity. Solid 
line shows linear extrapolation of the results for the conserved current. 
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FIG. 10. /tt in lattice unit as a function of rufa at M = 1.819 and Ng = 10 together with the 
experimental value. Data points for ZaA^^^ (filled squares) are slightly displaced horizontally for 
clarity. Solid line shows linear extrapolation of the results for the conserved current. 
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FIG. 11. /tt in lattice unit extrapolated to m/a 
with the experimental value. 



as a function of M aX Ng = 10 together 
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TABLES 



T'A RT T? T 


Siniulation parameters 


Oil io X oz X i\ g lattices 


in the quenched QCD. 


Ns 


4 


6 


10 


No. conf. 


20 


10 


30(M = 1.819), 20 


M 


1.7 


1.819 


1.7 




1.75 




1.819 




1.819 




1.9 




1.85 








1.9 








1.95 








2.0 






iTifa 


0.075 


0.075 


0.075 




0.050 


0.050 


0.050 




0.025 


0.025 


0.025 



-0.08 (M = 1.819) 

-0.086 (M = 1.819) 
-0.093 (M = 1.819) 
-0.10 (M = 1.819) 
-0.12 (M = 1.819) 



18 



